Impact of reionization on CMB polarization tests of slow-roll inflation 
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Estimates of inflationary parameters from the CMB B-mode polarization spectrum on the largest 
scales depend on knowledge of the reionization history, especially at low tensor-to-scalar ratio. 
Assuming an incorrect reionization history in the analysis of such polarization data can strongly bias 
the inflationary parameters. One consequence is that the single-field slow-roll consistency relation 
between the tensor-to-scalar ratio and tensor tilt might be excluded with high significance even if 
this relation holds in reality. We explain the origin of the bias and present case studies with various 
tensor amplitudes and noise characteristics. A more model-independent approach can account for 
uncertainties about reionization, and we show that parametrizing the reionization history by a set 
of its principal components with respect to E-moAe polarization removes the bias in inflationary 
parameter measurement with little degradation in precision. 



I. INTRODUCTION 

Temperature and polarization power spectra of the cos- 
mic microwave background (CMB) are consistent with 
predictions of the simplest inflationary models 
a nearly flat geometry, superhorizon correlations as 
probed by the spectrum of acoustic peaks, and primordial 
scalar perturbations that are adiabatic, Gaussian, and 
close to scale-invariant [1, H, 0, [1, • One of the key re- 
maining signatures of inflation, tensor perturbations (i.e. 
gravitational waves) [13, [Hi , has yet to be detected. De- 
pending on the amplitude of the tensor perturbations, 
which is not well constrained theoretically, it may be pos- 
sible to measure the angular power spectrum of the in- 
flationary gravitational waves in the i?-mode component 
of the CMB polarization on large scales. Non-detection 
of the tensor spectrum does not necessarily rule out in- 
flation, but upper limits on r can be used to exclude 
particular models of inflation and limit its energy scale. 
Many experiments have been proposed to search for this 

signal [12, [H [13, [H, [H, [13, m, [li, [13, [H • 

Measurement of tensor perturbations in the _B-mode 
polarization power spectrum would test models of infla- 
tion by constraining inflationary parameters. These pa- 
rameters include the tensor-to-scalar ratio, r, and the 
tensor spectral index, rit, which are related by a con- 
sistency relation under the simplest single-field slow-roll 
inflationary scenarios. If the tensor spectrum can be de- 
tected, precise measurements over a wide range of scales 
could test the consistency relation. 

CMB constraints on r and nt depend on the abil- 
ity to accurately determine the large-scale power in B- 
modes due to tensor perturbations, independent of the 
effects of other cosmological parameters. On the largest 
scales, the tensor B-mode spectrum depends not only 
on inflationary parameters but also on the reioniza- 
tion history of the universe [2^ . The main impact of 
reionization on the spectrum is through the total opti- 
cal depth, r. The 3-year WMAP measurements of E- 
mode polarization determine r to an accuracy of about 



30% I, [H, and future CMB experiments should con- 
strain r at the 5 - 10% level [H, [53, IH, [H] . 

However, the detailed evolution of the reionization his- 
tory also significantly affects the large-scale polarization 
spectra. Uncertainty in this history leads to added un- 
certainty in inflationary parameters. Moreover, incorrect 
inferences due to an oversimplified treatment of reion- 
ization may bias estimates of inflationary parameters. 
For unbiased estimation of the optical depth from the E- 
mode reionization peak, the solution is to use a complete, 
principal-component-based description of reionization ef- 
fects when estimating parameters from CMB polarization 
data [1^, [13]. In this paper, we extend this approach to 
tensor i?-mode polarization and show that it is equally 
if not more effective in ensuring accurate measurements 
with little loss in precision. 

The outline of the paper is as follows. We discuss the 
effects of reionization and inflationary parameters on the 
polarization power spectra and the large-scale degener- 
acy between these parameters in § |TT1 A brief overview 
of the principal component parametrization of the reion- 
ization history follows in § IIIII In § IIVI we describe our 
Markov Chain Monte Carlo analysis of simulated polar- 
ization data and give the resulting constraints on r, r, 
and rif, which we discuss further in §|Vl 



II. POLARIZATION DEPENDENCE ON 
REIONIZATION AND INFLATIONARY 
PARAMETERS 

Like the scalar _E-mode polarization power spectrum, 
the tensor i?-modc spectrum Cf^'"^ consists of two main 
components: one from recombination that peaks at € ~ 
100, and the other from the epoch of reionization, which 
peaks near ^ ^ 5 and dominates at £ < 20 (Fig. [T]). In 
the tensor power spectrum, these components arise from 
wavenumbers k ~ 0.01 Mpc"^ and 0.0007 Mpc"^ re- 
spectively. They provide a lever arm of over a decade in 
physical scale for measuring the tensor tilt. This lever 
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FIG. 1: _B-mode tensor spectra illustrating the degeneracy 
between r and nt for large-scale measurements, with angu- 
lar power spectra plotted in the upper panel and fractional 
deviations from the base model in the lower panel. For the 
base model (solid), r = 0.03, t = 0.1, and nt = —0.00375. 
The other two models have {r, nt} = {0.12, —0.00375} [long 
dashed) and {0.12,0.13} [short dashed), with a pivot scale 
of fcpivot ~ 0.01 Mpc~^. The reionization history here is as- 
sumed to be instantaneous. Cosmic variance of Cf ^'^ for 
the base model, which excludes the variance from lensing, is 
shown by the shaded band in the lower panel. 



arm is slightly shortened compared with the ratio of an- 
gular scales due to the closer distance to reionization. 
Nonetheless, these features provide an opportunity to 
test the slow-roll inflationary consistency relation [e.g.. 
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(1) 



between the tensor-to-scalar ratio r and tensor tilt nt- 
Deviations from the consistency relation and running of 
the tilt come in at second order in the slow-roll param- 
eters. Specifically, r is four times the ratio of the ten- 
sor power spectrum amplitude (for one component of 
gravitational waves) to the scalar curvature power spec- 
trum amplitude, in accordance with the definition used 
by CAMB and WMAP. 

Note that we quote r at a normalization scale of 
fcpivot = 0.01 Mpc~^ whereas 3-year WMAP results 
quote it at fcpivot = 0.002 Mpc^^ and CAMB defaults 
to fcpivot = 0.05 Mpc~^. Our choice, corresponding to 
the recombination peak at £ w 100, better refiects the 
scale at which the tensor spectrum can be best mea- 
sured. We will use subscripts to denote choices of scale 



other than 0.01 Mpc" 



'•g- ^0.05 foi' normalization at 



fcpivot = 0.05 Mpc , so that 



^fcpivot /Mpc 
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assuming no running of the scalar or tensor spectral in- 
dices. Since we only consider small deviations from scale- 
invariance, the tensor-to-scalar ratio varies little with the 
normalization scale so any corrections to Eq. ([1]) due to 
changing fcpivot are correspondingly small. 

Cosmic variance and reionization history uncertainties 
in the interpretation of the reionization peak limit the 
ability to measure the spectrum from this technique. In 
Fig. [l] we show that variation in the optical depth to 
reionization, r, can mimic changes to nt within the pre- 
cision of cosmic variance of the individual £ modes, 
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The example of Fig. [T] shows that there are at least 
two ways of breaking this degeneracy. One is by mea- 
suring the tensor spectrum at i > 100, where the tilt of 
the spectrum matters but the optical depth has no ef- 
fect other than an overall rescaling of the spectrum by 
e"^"^, which we absorb by changing the scalar amplitude 
As to keep AsC"^"^ fixed. Just as in the scalar case, if 
the tensor spectrum can be precisely measured beyond 
the recombination peak, the global constraint on nt will 
not be sensitive to reionization. Unlike the scalar case, 
Cf^'^ falls sharply at higher I and becomes masked by 
_B-modes generated by lensing of i?- modes [2^ [Cf^'^; 
see Fig. [2]). The power in the lensing S-modes is expected 
to be greatest a.t £ ^ 1000, a smaller scale than the reion- 
ization and recombination peaks of the tensor spectrum. 
The relative amplitude of the tensor and lensing contri- 
butions to Cf^ depends on the tensor-to-scalar ratio, 
which the 3-year WMAP data restrict to r < 0.3 

If r is near the current upper limits, the lensing spec- 
trum can be statistically subtracted to a large extent. 
Nevertheless, instrumental and foreground limitations 
may still prevent the extraction of information from 
scales beyond the recombination peak for next genera- 
tion experiments such as Planck jsol . [3ll . [s^ . [s^ . At much 
lower r, the best CMB constraints on nt will come from 
the combination of the recombination and reionization 
peaks unless the lensing signal can be subtracted directly 
from the polarization maps [13, [H, [s^, [13, [11] . 

The T — nt degeneracy is also broken through the con- 
straint on r from the E-mode reionization peak, which 
depends on r but not nt as long as r is small enough that 
Cf^ is dominated by scalar perturbations. To break the 
degeneracy in this way, it is important to have accurate 
constraints on r that do not depend on overly simplistic 
assumptions about the reionization history. 

The shape of the reionization peak depends on the 
history of the spatially-averaged ionized fraction, Xe{z). 
In Fig. [21 we illustrate two models with similar opti- 
cal depths but different reionization histories. A model- 
dependent analysis of the polarization power spectrum 
that uses an incorrect form of Xe{z) can result in sig- 
nificant bias in the total optical depth to reioniza- 
tion [H, [H, [H, [H. To the extent that the constraint 
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FIG. 2: E- and B-mode polarization angular power spec- 
tra and reionization histories (inset) for an extended, "dou- 
ble" reionization history with r = 0.090 (thick curves) and 
a nearly-instantaneous reionization model with optical depth 
r = 0.105 (thin). For the _B-modes, both the tensor spec- 
tra [with tensor-to-scalar ratio r = 0.3 (solid) and r = 0.03 
(dashed)] and lensing spectra are plotted. Dotted curves are 
the assumed Planck and CMBPol noise power spectra, from 
top to bottom at low I respectively. 

on fit relies on measurements of the reionization peak 
of Cf^'^, it may be biased as well in such an analysis. 
We shall see that the use of model-independent principal 
components of the reionization history protects against 
such biases at little cost to the precision of the inflation- 
ary test. 

III. PRINCIPAL COMPONENT 
PARAMETRIZATION OF REIONIZATION 

Following (25l . [27j , we parametrize the reionization his- 
tory as a free function of redshift by decomposing Xe(z) 
into its principal components with respect to the E'-mode 
polarization of the CMB: 

a;e(z)=a;f(z) + ^m^5^(z), (4) 

where the principal components, S'^(z), are the eigen- 
functions of the Fisher matrix that describes the depen- 
dence of Cf^ on Xe(z), are the amplitudes of the 
principal components for a particular reionization his- 
tory, and x^'^(z) is the fiducial model at which the Fisher 
matrix is computed. In practice, we construct the prin- 
cipal components assuming .E-mode information only by 
taking the tensor-to-scalar ratio r = 0. For allowed val- 
ues of r, the information on the ionization history from 
the tensors is subdominant to the scalars. The inverses 



FIG. 3: The five lowest-variance principal components of 
Xe(z) over redshifts & < z < 30, with increasing variance 
from top to bottom. 



of the eigenvalues of the Fisher matrix give the estimated 
variances of the principal components, cr^, which deter- 
mine the ordering of the components by requiring that 
cr^ < cr^_|_i. The main advantage of using principal com- 
ponents as a basis for Xe(z) is that only a small number of 
the components are required to completely describe the 
effects of reionization on large-scale CMB polarization, 
so we obtain a very general parametrization of the reion- 
ization history at the expense of only a few additional 
parameters. 

The principal components are defined over a limited 
range in redshift, Zmin < z < Zmax, with = at z > 

Zraax and = 1 at Z < Zmin- We take Zmin = 6, Zma x = 

30, and constant x^'^(z) = 0.15 here, and refer to [l^ 
for further discussion of the choices of these and other 
parameters related to the principal components. 

For complete representation of the effects of reioniza- 
tion between Zmin and Zmax on the low-^ E-mode spec- 
trum to better accuracy than cosmic variance, no more 
than the first five principal components are needed (as- 
suming Zniax ^ 30) [1^,[13]- Due to projection effects [4l|, 
the accuracy to which the lowest-variance principal com- 
ponents reconstruct Cf^'"^ at low I is even better than 
for the scalar Cf^. In the MCMC analysis presented 
in the following section, we always use the five lowest- 
variance principal components of Xe(z) with z^ax = 30, 
which we show in Fig. [3] (see [2^ for the effects of us- 
ing a different number of components to analyze E-mode 
data) . The amplitudes of these components then serve to 
parametrize general reionization histories in the analysis 
of CMB polarization data. 
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IV. MARKOV CHAIN MONTE CARLO 
CONSTRAINTS 



We use Markov Chain Monte Carlo (MCMC) analy- 
sis to explore the joint effects of the reionization his- 
tory and inflationary parameters on CMB polarization 
power spectra [see e.g.|42|,|4l,|43|. Chains of Monte Carlo 
samples are generated using the publicly available code 
CosmoMC m, which includes the code CAMB 
for computing theoretical angular power spectra at each 
point in the parameter space. We have modified both 
codes to allow specification of an arbitrary reionization 
history calculated from a set of principal component am- 
plitudes using Eq. as described in [25j . 

For chains in which x^iz) is parametrized by its prin- 
cipal components, the parameters that we vary include 
the amplitudes of the first five components, the tensor- 
to-scalar ratio at fcpivot = 0.05 Mpc^^, and the tilt of the 
tensor spectrum: {mi, TO2, 7713, TO4, 771,5, tq. 05, 'T'*}- When 
using the simple instantaneous form of Xe{z), the chain 
parameters are instead {t, ro.05, tij} with optical depth 
taking the place of principal components. We compute 
r at fcpivot — 0.01 Mpc""'^ as a derived parameter using 
Eq. ^ after generating each chain. 

When analyzing the parameter chains, we impose pri- 
ors on the principal component amplitudes corresponding 
to physical values of the ionized fraction. < Xe < 1 [25j . 
These priors are conservative in the sense that all ex- 
cluded models are unphysical, but the models we retain 
are not necessarily physical. 

We use only polarization for parameter constraints, 
and assume that the values of the standard ACDM pa- 
rameters (besides r) are fixed by measurements of the 
CMB temperature anisotropies. To a good approxima- 
tion the effect of Xe{z) on the large-scale polarization is 
independent of the other parameters [2^ . 

Specifically, we take f^^/i^ = 0.0222, n^h'^ ^ 0.106, 
1006* = 1.04 (corresponding to h ^ 0.73), A,e"^^ = 
1.7 X 10"^, and Us = 0.96, consistent with the WMAP 3- 
year temperature power spectrum. When computing the 
optical depth to reionization we take the helium frac- 
tion to be 1^ = 0.24 and assume that helium is neu- 
tral. With these parameter values, the reionization opti- 
cal depth out to Zmin = 6 is fixed at r(zinin) ~ 0.04. The 
remaining contribution to the total optical depth from 
Zmia < z < ^max IS determined by the values of {m^} for 
each sample in the chains. The default bin width for the 
fiducial models and principal components is Az = 0.25, 
which is small enough that numerical effects related to 
binning should be negligible. 

For each scenario that we study, we run four separate 
chains until the Gelman and Rubin convergence statistic 
R, corresponding to the ratio of the variance of param- 
eters between chains to the variance within each chain, 
satisfies R - 1 < 0.01 [13, [Hi. The convergence diag- 
nostic of [4^ is used to determine how much each chain 
must be thinned to obtain independent samples. 



A. Simulated data 

We use CAMB to generate model E- and i?-mode 
polarization power spectra. We take the data to be 
the exact values of a given model Cf^ = Cf^ and 
Cf^ ~ Cf^, neglecting cosmic variance and noise, so 
we expect constraints to be centered on the fiducial pa- 
rameter values rather than displaced by ~ 1 a. These 
constraints can be thought of as the average over many 
possible realizations of the data. (See [2^ for a discussion 
of the effects of cosmic variance when using realizations 
drawn from Ci instead of taking Ci as the data.) 

For the jth sample in a chain, the likelihood is 



lnL(,) = 



E 

1=2 



(5) 
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where Ci = Ct + Ng is the sum of the simulated data and 
noise spectra, and C^q) = C£(-,) -|- Ng is the theoretical 
spectrum calculated with the parameter values at step 
j in the chain plus the noise spectrum. For consistency 
with CosmoMC, the likelihood contains an extra factor 
of /sky , the fraction of sky observed [s^ . 

We set A'^ = to simulate measurements limited only 
by cosmic variance. For more realistic scenarios based on 
Planck and CMBPol, we model noise as 



.-1/2 



Nf 



/.jK-rad 

-1/2 



exp 



FWHM 



/rad)2 



;ln2 



(6) 



where Wp is the polarization noise level and ^fwhm is 
the beam width. 

We typically compute the likelihood up to ^max = 1000. 
The i?-mode spectrum and B-mode tensor spectrum are 
fixed at ^ ^ 100 by setting Ase~'^'^ constant in the Monte 
Carlo chains. The impact on r and rit constraints of the 
_B-mode spectrum at niultipoles greater than a few hun- 
dred is negligible in the presence of the lensing spectrum 
and experimental noise for all allowed values of r. 

The fiducial reionization history is the extended, dou- 
ble reionization model with polarization spectra plotted 
in Fig. [21 for which r = 0.090. This function is cho- 
sen because it is not well described by the instantaneous 
reionization model, so biases due to assuming instanta- 
neous Xe{z) should be readily apparent in the parameter 
constraints. 

We consider two fiducial values of the tensor-to-scalar 
ratio, To. 05 = 0.3 and ro.05 ~ 0.03. Using the consistency 
relation [Eq. ^] to set rit and taking the fiducial scalar 
tilt Hs = 0.96, these tensor-to-scalar ratios correspond 
to r = 0.299 and r = 0.0283 at fcpivot = 0.01 Mpc^^, 
respectively. 

Assuming a power-law primordial spectrum, the larger 
fiducial tensor-to-scalar ratio, 7'o.o5 ~ 0.3, is approxi- 
mately the 68% upper limit on r from 3-ycar WMAP 
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TABLE I: Constraints on r, r, and nt for simulated data based on the double reionization history with optical depth r — 0.09. 
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0.0288 ± 0.0006 
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0.03 


0.0283 


-0.00375 


CMBPol 


Yes 


0.092 ± 0.007 


0.032 ± 0.009 


0.069 ±0.184 


0.03 


0.0283 


-0.00375 


CMBPol 


No 


0.122 ± 0.005 


0.038 ±0.010 


0.491 ±0.165 



data, although when large-scale structure data are in- 
cluded it is ruled out at about 95% CL 0. A ten- 
sor spectrum with this value of r should be detectable 
by the Planck satellite [l^. The smaller value of r 
may be out of the reach of Planck but accessible to 
a next-generation CMB satellite, such as the proposed 
CMBPol In addition to CV-limited measurements, 
we consider a noise spectrum based on Planck with sen- 
sitivity Wp =81 ^iK' and beam size ^fwhm = 7.1' 
for the To. 05 — 0.3 simulated data (l^ . [50| . and low- 
resolution CMBPol-Wke noise with Wp =20 /iK' and 
^^FWHM = 60' for To. 05 = 0.03 [ll] (dotted curves in 
Fig. [2]). We assume /sky = 0.8 for both Planck and 
CMBPol and take /sky = 1 for the more idealized, CV- 
limited data. In all cases we neglect foregrounds, as- 
suming that they can be adequately subtracted using 
polarization data from frequency channels not used for 
cosmological parameter estimation. 

The i?-mode lensing spectrum has a significant impact 
on constraints on r and nt, but including lensing signifi- 
cantly slows down the computation of the angular power 
spectra in CAMB. However, Cf^'^ is nearly indepen- 
dent of the parameters that we vary in the Monte Carlo 
chains. Rather than computing the effects of lensing di- 
rectly, then, we treat the lensing spectrum as a fixed con- 
tribution to the noise power spectrum. Tests comparing 
this approximation to an analysis including the full ef- 
fects of lensing show that the constraints obtained for r 
and nt are the same. Using r-dependent lensing spec- 
tra computed for each Monte Carlo sample appears to 
improve the constraint on r, but this is an artifact due 
to fixing AsB"^"^ and other parameters that also affect 
Cf^'^, mainly 0,„/i2 ^ 

B. Results 

Table U lists the ID marginalized constraints on r, r, 
and nt from the MCMC analysis for each case study. 
The constraints on r, r, and nt are plotted in Figs. [3] 
and [H] as 2D contours after marginalizing over all other 



parameters. The two sets of contours in each panel use 
the same simulated data, but different parametrizations 
of Xe{z). For the thick contours, the five lowest- variance 
principal components of Xe{z) are included in the Monte 
Carlo chains, while the thin contours come from chains 
that treat Xe{z) as instantaneous reionization with only 
one parameter, r. For the principal component chains, 
T is derived from the principal component amplitudes, 

Since the constraint on optical depth comes primarily 
from the reionization peak of the _E-mode spectrum, es- 
timates of r in all cases in Tab. |T] are affected similarly by 
the parametrization of the reionization history. As noted 
in [3^ [5^. using the one-parameter instantaneous Xe(z) 
when that model is not sufficient to describe the true 
reionization history (in this case, the double reionization 
history represented by the thick curves in Fig.[2|) can lead 
to a significant bias in the constraint on r. Using a set of 
principal components to parametrize Xe{z) removes the 
bias in the optical depth f2^. 

The impact of the reionization history assumptions for 
infiationary parameters depends more strongly on the 
values of r and the noise. We consider first the most 
optimistic scenario for measurement of the _B-mode ten- 
sor spectrum: a CV-limitcd measurement of Cf^'^ with 
approximately the largest amplitude currently allowed, 
''0.05 = 0.3. With this tensor-to-scalar ratio, the tensor 
spectrum dominates over the B-mode lensing spectrum 
for £ % 150 (see Fig. [2|), so the recombination peak of 
^BB,T fj-Qjjj contamination over roughly a decade 

in £. Moreover, the lensing spectrum may be statisti- 
cally subtracted well beyond the multipole where it con- 
tributes equal power. 

Since the amplitude of the tensor spectrum depends on 
T, r, and nt, we expect some degeneracy between these 
parameters, as described in § |TT1 Because of this degen- 
eracy, a bias in r can generate biases in the infiationary 
parameters as well, mainly nt (since r at fc = 0.01 Mpc~^ 
is tightly constrained by measuring Cf""'' at ^ « 100). 
However, r only affects the \ow-£ reionization peak (with 
our assumption that Ase~'^'^ is held fixed), so these biases 
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FIG. 4: 2D marginalized 68 and 95% contours for the opti- 
cal depth (r), tensor-to-scalar ratio (r), and tensor spectral 
index (nt). For the simulated polarization spectra, we take 
Xe{z) to be a double reionization history with r — 0.090 (see 
Fig- nil- The fiducial tensor-to-scalar ratio is set to ro.os ~ 0.3 
(r = 0.299 at A; = 0.01 Mpc"^), and the fiducial tensor tilt is 
assumed to obey the consistency relation, nt = — r/8 (right 
panels, line). Crosses indicate these fiducial parameter val- 
ues. For the thick contours, the chains include the five lowest- 
variance principal components of Xe{z). Thin contours show 
the constraints when Xe{z) is instead assumed to be instanta- 
neous and parametrized only by r. The lensing spectrum is 
treated as a contribution to the noise. Top panels: CV-hmited 
data; bottom: including Planck-like noise. 



will only be significant if the reionization peak con- 
tributes significantly to constraints on the tensor spec- 
trum. In this optimistic scenario with tq.os = 0.3 and cos- 
mic variance-limited measurements, enough of the ten- 
sor recombination peak is observable so that additionally 
measuring the reionization peak has a negligible effect on 
constraints on r and nt; the values of these parameters 
are determined almost entirely by Cf^ at 20 ^ ^ < 500. 
As a result, the bias in the infiationary parameters due to 
incorrectly assuming instantaneous reionization is small; 
as the contours in the top right panel of Fig. [4] show, the 
true parameter values remain within the 68% CL even 
with this bias. 

Now consider the same tensor-to-scalar ratio but more 
realistic assumptions about the experimental noise. Us- 
ing our assumed noise spectrum for Planck as described 
in § IIV Al we find constraints on t, r, and nt as shown 
in the bottom panels of Fig. U) The most obvious dif- 
ference from the CV-limited case is that the parameter 
uncertainties are larger, especially for the inflationary pa- 



FIG. 5: Same as Fig.lH but with fiducial tensor-to-scalar ratio 
ro.05 = 0.03 (r = 0.0283 at A; = 0.01 Mpc"^). Top panels: 
CV-limited data; bottom: including CMBPol-like noise. 



rametcrs. There is also a difference in the instantaneous 
reionization biases: while the optical depth bias is simi- 
lar to the CV-limited case, nt is biased much more with 
Planck-like noise than with a CV-limited measurement. 
The additional bias is due to the greater dependence on 
the reionization peak of Cf^'"^ for nt constraints since 
the noise associated with Planck makes measurement of 
the recombination peak much more difficult. 

Note that while r at fcpivot = 0.01 Mpc~^ is still de- 
termined fairly accurately in this scenario, if the tensor- 
to-scalar ratio were quoted elsewhere, e.g. as ro.oo2 or 
^0.05 1 £^ significant bias would appear in that parameter 
as well. For example. Fig. [5] shows that ro.oooTj normal- 
ized at the approximate scale of the reionization peak, 
lies significantly below the fiducial value in most of our 
test cases. The choice of scale also affects the uncertainty 
in the tcnsor-to-scalar ratio; while the error we quote for 
ut agrees with similar studies, e.g. [sij . the error on r dif- 
fers in general since degeneracy with nt due to a different 
choice of fcpivot can greatly increase Ur ■ 

With the larger bias in nt in the presence of Planck- 
like noise, the true parameter values are excluded at the 
95% CL. Moreover, the consistency relation for single- 
field slow-roll inflation (plotted as a line in the r — nt 
plots) is excluded at about the same confidence level. An 
incorrect assumption about the reionization history could 
therefore lead to wrongly rejecting the simplest class of 
inflationary models. 

As the thick contours in Fig. 3] show, when general 
reionization histories are considered by parametrizing 
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0.1 0.2 0.3 



0.01 0.02 0.03 



FIG. 6: ID marginalized distributions of ro.0007, the tensor- 
to-scalar ratio on scales near the reionization peak, for each 
of the four MCMC scenarios in Figs. |4] and [5] The fiducial 
values of ro.oooT are indicated by vertical dashed lines, and 
the normalization is arbitrary. 



Xe{z) by its principal components, the biases in r, r, and 
nt are all removed. The constraints on r and rit still may 
not be very strong for an experiment like Planck, but at 
least when using the principal components of Xe{z) one 
would not be led to exclude the true underlying model of 
inflation. Furthermore, the precision of the constraints 
is not substantially degraded by the principal component 
parametrization. 

Note that the exact magnitude and direction of the 
parameter biases depend on what the true reionization 
history actually is. The goal here is not to make specific 
predictions about these biases, but rather to give an idea 
of how large they could be and to show that they can be 
avoided by allowing for general models of reionization. 

How does uncertainty about reionization influence con- 
straints on inflationary parameters if the amplitude of the 
tensor spectrum is smaller? In general, we expect the 
importance of the reionization peak of Cf^""^ relative to 
the recombination peak to increase as the tensor spec- 
trum drops further below the dominant source of "noise" 
on large scales, whether that be the noise spectrum of 
an experiment or the i?-mode lensing signal. To check 
this expectation, we use simulated data with a fiducial 
tensor-to-scalar ratio of ro.05 = 0.03. The resulting con- 
straints on r, r, and nt, plotted as for the ro.05 = 0.3 
results, arc shown in Fig. [5] for both cosmic variance- 
limited measurements (top panels) and the CMi?Po/-like 
noise spectrum described in IIV Al (bottom panels). 

With a CV-limited experiment and ro.05 = 0.03, due 
to the greater influence of lensing relative to the tensor 
spectrum, the reionization peak < 20) and recombina- 
tion peak at ^ ^ 300 contribute roughly equally to the 



constraint on rit. Because of this, assuming the wrong 
reionization history biases nt more than for ro.05 = 0.3, 
and in this case the consistency relation would be ex- 
cluded at more than the 95% CL. 



If we add the CMBPol noise spectrum {wp = 
20 /iK', 0FWHM — 60'), both biases and errors on pa- 
rameters are larger. The constraints in this case are 
quite similar to the constraints in the case of Planck with 
''0.05 = 0.3 (except with all r values reduced by a factor 
of 10), which makes sense since the low-£ amplitude of 
CMBPol noise relative to ' with ro.05 = 0.03 is sim- 
ilar to that of Planck noise to Cf^'"^ with ro.05 = 0.3 (see 
Fig. [21). Note that our assumptions about the CMBPol 
noise spectrum are on the conservative end of the range 
usually considered for such an experiment; better sensi- 
tivity (~ 1 /iK') and/or resolution 1 — 10 arcmin) 
would enable measurement of the tensor i?-modcs to 
smaller scales so that constraints on inflationary param- 
eters would be closer to the CV-limited scenario. 



For even lower values of r than those considered here, 
the impact of reionization on constraints on inflationary 
parameters is likely to be the same or greater. Even 
for an idealized, cosmic variance-limited experiment, the 
B-mode signal due to lensing becomes a significant con- 
taminant of Cf^'"^ at r < 0.03. Due to the shape of the 

spectra, if any part of Cf^'^ is detectable above the lens- 
ing S-modes and noise for low tensor-to-scalar ratios it 
will be the reionization peak. Any bias in r due to incor- 
rect modeling of the reionization history will then cause 
the inferred value of nt to be biased. Of course, at very 
low r the uncertainties are so large that the parameters 
are not usefully constrained at all. In this case, it may be 
possible to measure nt using direct observations by future 
gravitational-wave experiments to better accuracy than 
what is possible with CMB polarization [HI, [HI] . 



Finally, note that although the choice of normaliza- 



tion scale fcpivot = 0.01 Mpc^^ is intended to decorrelate 
r from the other parameters, there is some remaining 
degeneracy with nt- The direction of this degeneracy 
depends on whether or not Cf^''^ can be measured on 
scales smaller than fcpivot- For CV-limited measurement 
of a i?-mode spectrum with r = 0.3, these smaller scales 
are observable. A larger tensor tilt that increases the 
power at ^ ^ 100 can be compensated for by lowering 
the value of r, so r and nt arc anticorrelated in this case 
(top right panel of Fig. |4|). If the tensor spectrum on 
such small scales is hidden by lensing B-modes or noise, 
then the low-^ side of the recombination peak becomes 
more important. Increasing tensor tilt lowers the power 
at 20 < ^ < 100, so r is correlated with nt so as to 
match the spectrum of the data on these scales (lower 
right panel of Fig. 2] and right panels of Fig. [5]). 
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V. DISCUSSION 

The value of the optical depth to reionization esti- 
mated from the CMB i?-modc polarization spectrum on 
large scales can be biased by adopting a model that has 
insufficient freedom to describe the true reionization his- 
tory. Likewise, the use of simple reionization models can 
bias inflationary parameters such as the tensor-to-scalar 
ratio and tensor tilt that depend on the large-scale ampli- 
tude of the i3-mode spectrum of primordial gravitational 
waves. In each case, the problem can be solved by using a 
more general parametrization of the reionization history. 
We have shown that using a small but complete set of 
the principal components of the reionization history ef- 
fectively yields unbiased constraints on both reionization 
and inflationary parameters. 

Measurements of r and rit are only affected by the as- 
sumed form of the reionization history if the reionization 
peak of the tensor i3-mode spectrum at the very largest 
scales is needed to precisely constrain the parameters. If, 
instead, good constraints can be obtained using only the 
_B-mode recombination peak at intermediate scales, then 
assumptions about reionization do not affect tests of the 
consistency relation between r and rit- They would in- 
stead appear as false evidence for running of the tensor 
tilt in violation of slow-roll expectations. Measurement 
of the recombination peak however is inhibited by experi- 
mental noise and contamination from i?-mode power con- 
verted to i?-mode power by gravitational Icnsing, both of 
which become more important at smaller scales. 

To study the potential impact of reionization on pa- 
rameter constraints from _B-mode polarization, we have 
employed a Markov Chain Monte Carlo analysis of sim- 
ulated CMB polarization power spectra and compared 
results for two descriptions of reionization: a simple, 
one-parameter, instantaneous reionization model, and a 
parametrization using principal components of the reion- 
ization history with respect to the E-mode polarization 
power spectrum. 

By varying the properties of the simulated polariza- 



tion power spectra, including the fiducial tensor-to-scalar 
ratio and the noise spectrum, we have determined over 
what range of scales CMB polarization data is most im- 
portant for constraining inflationary parameters in vari- 
ous scenarios. In particular, the question of whether the 
large-scale reionization peak of Cf^'"^ or the smaller- 
scale recombination peak is more important determines 
the severity of bias in inflationary parameters when reion- 
ization is modeled incorrectly. 

If the tensor-to-scalar ratio is near the current upper 
limit oi r ^ 0.3 and measurements of B-mode polariza- 
tion are limited only by cosmic variance, then the spec- 
trum on scales 20 ^ ^ ^ 500 dominates constraints on 
r and nt and incorrect assumptions about reionization 
do not strongly bias the results. If the true tensor-to- 
scalar ratio is more than a factor of a few smaller than 
this upper bound, however, then lensing _B-modcs limit 
the information that can be extracted from the recombi- 
nation peak of the tensor spectrum alone. Furthermore, 
all-sky experiments in the foreseeable future are likely to 
have noise that exceeds the lensing signal, making tests 
of inflation even more reliant on the reionization peak 
of the tensor B-modes on large scales. In all of these 
cases, a general parametrization of reionization such as 
that provided by principal components allows the use of 
the i?-mode reionization peak for inflationary parameter 
constraints without significantly worsening the errors on 
those parameters. 
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